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Fourth-order cumulants of physical quantities have been used to charac- 
terize the nature of a phase transition. In this paper we report some Monte 
Carlo simulations to illustrate the behavior of fourth-order cumulants of mag- 
netization and energy across second and first-order transitions in the phase 
diagram of a well known spin-1 Ising model. Simple ideas from the theory 
of thermodynamic fluctuations are used to account for the behavior of these 
cumulants. 

I. INTRODUCTION 

There are many attempts to characterize the order of a phase transition on the basis 
of the analysis of numerical data obtained from simulations of finite spin systems. One 
of the approaches to this problem consists in the analysis of the behavior of fourth-order 
cumulants of physical quantities (as the order parameter and the energy) associated with the 
systems under consideration Properties of the fourth-order cumulants of magnetization 
(and energy) have been investigated in the context of finite-size effects in magnetically (and 
thermally) driven first-order transitions in Ising and Potts models (as well as in the 
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case of some other systems 

In this paper, we perform Monte Carlo simulations for the well known Blume-Capel 
model to illustrate the behavior of the fourth-order cumulants of magnetization and 

energy across first and second-order transitions in the phase diagram of this system. We 
show that it is possible to draw some conclusions from the study of relatively small lattices. 
The general features of the cumulants can be accounted for by simple arguments from the 
theory of thermodynamic fluctuations. In particular, we emphasize the differences between 
the two types of cumulants, and the alternative definitions of the cumulant of energy (which 
has not been fully appreciated in previous investigations). 

The layout of this paper is as follows. In Section 2 we define the cumulants of a physical 
quantity. In Section 3 we introduce some further definitions, and discuss some properties of 
the Blume-Capel model. Simulations for the fourth-order cumulants of magnetization and 
energy are reported in Sections 4 and 5, respectively. We hope to have provided another 
example of the use of these cumulants to characterize the order of a phase transition. 

II. DEFINITION OF THE CUMULANTS 

The cumulants of a quantity x can be obtained from an expansion of the form 

1 1.1 

< exp (x) >= 1+ < x > H — < x 2 > H — < x 3 > H < x 4 > +... , (1) 

2 6 24 



where the brackets denote an average [TT| • If we keep terms up to fourth order, the logarithm 
of this expansion may be written as 

ln{< exp (z) >} =< x > +ig 2 + ^Q 3 + ^Q 4 + - , (2) 

where the cumulants Q 2 , Q3, and Q 4 , are given by 

Q 2 =< x 2 > - < x > 2 , (3) 

Q 3 =< x 3 > -3 < x >< x 2 > +2 < x > 3 , (4) 
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and 

g 4 =< a; 4 > -3 < x 2 > 2 -4 < x >< x 3 > +12 < x > 2 < x 2 > -6 < x > 4 . (5) 
These cumulants can be rewritten in the more compact form, 

Q 2 =< (x- < x >) 2 >, (6) 

Q 3 =< (x- < x > f >, (7) 



and 



Q 4 =< (a;_ < x >) 4 > -3 < (x- < x >y > 2 . (8) 



Also, the fourth-order cumulant Q± is more often written as 



< (x- < x > L ) A > 



v -w - 1 - 3 ;\;-;;v 1 v>v (9) 

where L is the linear size of the lattice under consideration. 

III. THE BLUME-CAPEL MODEL 

The Blume-Capel model is given by the spin Hamiltonian 

N N 

n = - J E s i s i + Dj2s 2 -Hj2s t , (io) 

i=\ i=l 

where Si = +1, 0, —1, on sites i = 1, N of a Bravais lattice, and the first sum is performed 
over all pairs of nearest-neighbor sites. We consider the ferromagnetic case, with positive 
exchange (J > 0) and anisotropy (D > 0) parameters, which gives rise to a competition 
between distinct spin orderings. In the D/J versus T/J space, where T is the absolute 
temperature, the phase diagram consists of an ordered and a disordered phase separated by 
a transition line that changes character from first to second-order at a well defined tricritical 
point. We use this model to illustrate the behavior of the fourth-order cumulants. 
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In zero field, the ordered phase is characterized by symmetric magnetizations, +m and 
— to , with the same energy. In general, we have (H) L ^ 0, and even (H n ) L ^ 0, for all n. 
Therefore, the fourth-order cumulant of energy is written as 

Ve{l) ~ 1 " 3< (H-<H> L y >r ( } 

where L is the linear size of the lattice. As the magnetization is symmetric, that is, < 
to™ >l= for n odd, the fourth-order cumulant of the magnetization is given by 

It should be remarked that, as < E n for all n (including odd values of n), and 

< to™ >l= for odd n, the energy and the magnetization give rise to distinct expressions 
of the fourth-order cumulant. 

Now we perform a preliminary Monte Carlo simulation to look at the form of the distri- 
bution of probabilities for the magnetization (and to motivate the choice of the distributions 
of energy and magnetization to be used in the forthcoming theoretical calculations). Let us 
define the dimensionless variables 

T , D , , H 

t = — , a = — , and h = — , (13) 

and use a Metropolis algorithm to perform Monte Carlo simulations for the Blume-Capel 
model, in zero field, on a simple cubic lattice of side L. 
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t = 0.1667 , lattice 8x8x8 
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FIG. 1. Histogram of the absolute values of the magnetization, \m\, across a first-order phase 
transition. 

In Fig. (1), we show the distribution of the absolute values of the magnetization for 
t = 0.1667, L = 8, and several values of d. In the top left graph, for d = 0.5000, \m\ has a 
unique maximum around \m\ = 0.9, which corresponds to an ordered phase. The bottom 
left graph, for d = 0.5008, displays two peaks (around m = and \m\ = 0.9), which indicate 
the coexistence of phases (m = 0, and m = ±m ^ 0). Upon increasing the value of d, the 
peak at \m\ = is enhanced, while the peak around \m\ = 0.9 is depressed (see the top 
right graph, for d = 0.5017). For even larger values of d, there remains a single peak around 
\m\ = 0, which indicates that the system is in the paramagnetic phase (see the graph for 
d = 0.5050). This set of graphs represents a first-order transition. 
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d = 0.35 , lattice: 8x8x8 

0.10 



t = 0.3667 




0.8 1 




0.6 0.8 10 



t = 0.3833 




0.2 0.4 0.6 0.8 1 

Iml 



0.10 



1=0.06 




t = 0.5000 



0.2 0.4 0.6 0.8 1.0 

Iml 



FIG. 2. Histogram of the absolute values of the magnetization, \m\, across a second-order phase 
transition. 

In Fig. (2), we show a histogram of \m\ across a second-order transition, for d = 0.35. 
The top left graph, for t = 0.3667, corresponds to the ordered phase, with a single peak of 
the absolute value of the magnetization at |m| = 0.5. As the value of t is increased, this peak 
moves toward \m\ = and no other peak arises (see the lower left graph, for t = 0.3833). 
For higher values of t, \m\ peaks at \m\ = (see the graphs at right, for t = 0.4167 and 
t = 0.5000). This set of graphs illustrates a continuous phase transition. 

Figs. (1) and (2) provide the motivation for choosing a Gaussian form for the probabil- 
ity of magnetization, p(m). Although we are showing data for the magnetization, similar 
histograms can be built for the energy, which also give support to a Gaussian-shaped prob- 
ability distribution, p(E). For large lattices [L — > oo) these distributions are expected to 
tend to Dirac delta functions. The thermodynamic consistency of these assumptions has 
been discussed in detail by Challa, Landau, and Binder J| (see also the work of Oitmaa and 
Fernandez ||12||). 
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Now we present separate analyses of the fourth-order cumulants associated with magne- 
tization and energy across first and second-order phase transitions. 



IV. CUMULANTS OF MAGNETIZATION 



In the ordered phase the distribution of probability of the magnetization consists of two 
peaks around +m and — m . For small lattices (L finite), we assume the Gaussian form 



p(m) = -C exp 
2 



(m + m o y 
2^ 



-C exp 



2a 2 



(14) 



where C = (27rcr 2 ) _1//2 is a normalization constant, and the parameter a should be inversely 
proportional to the lattice volume, L 3 . In the L — > oo limit, we have two Dirac delta 
functions, 



p(m) — > ^5(m + m ) + ^5(m - m ) 



(15) 



For finite lattices, with p(m) given by Eq. ([14]) , we have < m n >= 0, for odd values of 



n, and 



< m >- 



m 2 p{m)dm = a 2 + m 2 , 



(16) 



and 



< m >- 



m p{m)dm = 3a + 6cr m D + m G 



(17) 



Thus, the fourth-order cumulant is given by 



3a A + Qa 2 m 2 o + m 4 



(18) 



In the disordered phase, p(m) has a single peak at m = 0. Inserting mo = in Eq. ([T8"D , 
we have V m (L) = for the disordered phase. 

For infinite lattices (a — > 0), it is easy to see that V m (L) — > 2/3 in the ordered phase. In 
the disordered phase, V^(L) — > 0, regardless of the value of the parameter a. These limiting 
values for the infinite lattice also come from the double-delta distribution p(m), given by Eq. 



(|15l). In this case, in the ordered phase, we have < m n >= 0, for odd n, and < m n >= m™ 
for even n. Thus 
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FIG. 3. Fourth-order cumulant of the magnetization versus temperature, tq, across a sec- 
ond-order phase transition. Monte Carlo data were obtained for the Blume-Capel model on a 
cubic lattice (q = 6) of side L = 8, for d = 0.15. Averages were calculated from 50000 Monte Carlo 
steps after thermalization. 

In a second order phase transition, the two peaks of the distribution of probabilities 
p(m) in the ordered phase move towards each other and form a unique peak at m = as the 
system passes to the disordered phase. Fig. (3) shows the cumulant of magnetization for the 
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Blume-Capel model in a second order phase transition. These simulations were performed 
for a cubic lattice of side L = 8 (and coordination q — 6), at d — 0.15. We have run 5000 
times through the lattice to reach thermalization. Each average was then calculated using 
50000 additional steps. 

In a first-order phase transition we have the coexistence of the ordered and disordered 
phases. The distribution of probabilities p(m) has peaks at m = ±m ^ 0, and m = 0. For 
an infinite lattice, we take the triple-delta distribution, 

p{m) = c5{m + m ) + c5{m — m ) + (1 — 2c)<5(m) , (20) 

where c is a positive constant. We then have 

/oo 
m 2 p(m)dm = 2cm 2 , (21) 
-oo 

and 

/oo 
m A p(m)dm = 2cm 4 . (22) 
-oo 

Therefore, 

V -M = 1 -3®^ = 1 -k- (23) 

For small values of c, namely c < 1/6, the cumulant V m (L) is negative. Similar results could 
have been obtained from the L — > oo limit of a Gaussian form of p{m) for a finite lattice. In 
Fig. (4), we show the fourth-order cumulant of the magnetization across a first-order phase 
transition. The simulations were performed for a cubic lattice, with L — 8, at t — 0.15. 
We have run 5000 steps through the lattice to reach thermalization. Each average was then 
calculated from 50000 additional steps. In the ordered phase, we do have V m (L) =2/3. In 
the disordered phase there are very small fluctuations around V m (L) = 0. At the transition 
V m (L) assumes a negative value. 
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FIG. 4. Fourth-order cumulant of the magnetization versus the anisotropy, dq, across a first-order 
phase transition. Monte Carlo data were obtained for the Blume-Capel model on a cubic lattice 
(q = 6) of side L = 8, at temperature t = 0.15. Averages were calculated from 50000 Monte Carlo 
steps after thermalization. 



V. CUMULANTS OF ENERGY 

In both ordered and disordered phases, the probability distribution of energy has a unique 
peak at a certain value, which we call E Q . For small lattices, we can write the Gaussian form 

{E-E Q f- 



p(E) = Cexp 



2a 2 



(24) 



where C = (27ra 2 ) 1 ^ 2 . For large lattices (L — > oo), we have 



p(E)^S(E-E ). 



(25) 
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^From Eq. (|24|), we obtain 

/oo 
Ep(E)dE = E , (26) 
-co 

/oo 
(E - E ) 2 p(E)dE = a 2 , (27) 
-co 



and 



/CO 
(E - E fp(E)dE = 3cr 4 . (28) 
-co 



Inserting these expressions into Eq. (|TT|) , we have 



V ^L) = 1 - = 0, (29) 



for all values of the parameter a (that is, independently of the size of the lattice used in the 
simulation). Using the limiting distribution, given by Eq. (^), we also have 

/CO 
E5(E - E )dE = E , (30) 
-co 

and 

/CO 
(E - E ) n 6(E - E Q )dE = 0, (31) 
-co 

for all n. Hence, Ve(L) —>■ 0, as we have already obtained. 
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FIG. 5. Fourth-order cumulant of energy versus temperature, tq, across a second-order phase 
transition. Monte Carlo data were obtained for the Blume-Capel model on a cubic lattice (q = 6) 
of side L = 8, for d = 0.15. Averages were calculated from 50000 Monte Carlo steps after thermal- 
ization. 

In a second-order transition the distribution p(E) displays just a single peak, that moves 
from an initial value E\ to a final value E2. Therefore, the cumulant of energy across a 
second-order transition always vanishes, independently of the lattice size. The Monte Carlo 
estimates of Ve(L) for the Blume-Capel model, as shown in Fig. (5), for d — 0.15 and lattice 
size L = 8, indicate a small maximum next to a small minimum near the second-order phase 
transition. This behavior suggests that it becomes too simple to describe the probabilities 
in the immediate neighborhood of a continuous transition by a symmetric Gaussian form. 
This is also hinted by the bottom left graph of Fig. (2), which is already quite asymmetric. 
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In a first-order phase transition there is a coexistence between the ordered phase, associ- 
ated with a distribution of energy peaked at Ei, and the disordered phase, with a distribution 
peaked at E 2 . Then, we can write 

p(E) = c5(E-E 1 ) + (l-c)6(E-E 2 ), (32) 

from which we have 

/oo 
Ep(E)dE = cE 1 + (1 - c)E 2 , (33) 
-oo 

/oo 
[E - cE x - (1 - c)E 2 fp(E)dE = c(l - c){E 1 - E 2 ) 2 , (34) 
-oo 

and 

/oo 
[E - cE 1 - (1 - c)E 2 fp(E)dE = 
-oo 

= c(l - c) (1 - 3c + 3c 2 ) (^ - E 2 ) 4 . (35) 
The fourth-order cumulant is given by 

which becomes negative for small values of either c or 1 — c. In Fig. (6), for a cubic lattice 
of size L — 8, at t — 0.15, we show the fourth-order cumulant associated with the energy of 
the Blume-Capel across a first-order phase transition. 
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FIG. 6. Fourth-order cumulant of energy versus the parameter of anisotropy, dq, across a 
first-order phase transition. Monte Carlo data were obtained for the Blume-Capel model on a 
cubic lattice (q = 6) of side L = 8, at temperature t = 0.15. Averages were calculated from 50000 
Monte Carlo steps after thermalization. 

In Figs. (7) and (8), we illustrate the fourth-order cumulants of magnetization and energy 
across a second-order transition. We see that the cumulants V m (L), for different values of L, 
cross at a unique point, which can be used to estimate the transition temperature [p|,|6|,|?|,|T3| . 
However, it should be pointed out that a precise location of the transition requires a detailed 
study of finite size scaling, which is beyond the scope of this paper |Q-42,||,[I^-|T^] . 
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FIG. 7. Fourth-order cumulant of magnetization versus temperature, tq, across a second-order 
phase transition (d = 0.15), for the Blume-Capel model on a sequence of cubic lattices (q = 6), of 
lattice sizes L = 4, 8, and 14. The Monte Carlo averages were calculated from 50000 lattice steps 
after thermalization. 
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FIG. 8. Fourth-order cumulant of energy versus temperature, tq, across a second-order phase 
transition (d = 0.15), for the Blume-Capel model on a sequence of cubic lattices (q = 6), of lattice 
sizes L = 4, 8, and 14. The Monte Carlo averages were calculated from 50000 lattice steps after 
thermalization. 

Many authors use a fourth-order cumulant of the energy given by the form 



< E 4 >i 

' L 



instead of the connected expression of Eq. (|PT|). Although this may work for Ising and Potts 



models ||17|| , it is important to emphasize that for the Blume-Capel model we have to use the 
correct definition, given by Eq. (|HD , to be able to extract the order of the phase transition 
(see also the recent works of Janke [18|, and of Borgs and collaborators ||16|| ). 
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VI. CONCLUSIONS 



We have used the spin-1 Ising model of Blume and Capel to illustrate the feasibility 
of characterizing the order of a phase transition from a simple analysis of the behavior 
of the fourth-order cumulants of energy and magnetization. The general features of these 
cumulants can be derived from simple arguments of the theory of thermodynamic fluctua- 
tions. In the literature, there are two definitions of the fourth-order cumulant of energy. We 
have pointed out that, in the case of systems as the Blume-Capel model, it is important to 
consider the connected form of the cumulant of energy. 
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